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We present the optical spectra of Gei-^Si^ alloy nanocrystals calculated with time-dependent 
density-functional theory in the adiabatic local-density approximation (TDLDA). The spectra 
change smoothly as a function of the composition x. On the Ge side of the composition range, the 
lowest excitations at the absorption edge are almost pure Kohn-Sham independent-particle HOMO- 
LUMO transitions, while for higher Si contents strong mixing of transitions is found. Within TDLDA 
the first peak is slightly higher in energy than in earlier independent-particle calculations. However, 
the absorption onset and in particular its composition dependence is similar to independent-particle 
results. Moreover, classical depolarization effects are responsible for a very strong suppression of 
the absorption intensity. We show that they can be taken into account in a simpler way using 
Maxwell-Garnett classical effective-medium theory. Emission spectra are investigated by calculat- 
ing the absorption of excited nanocrystals at their relaxed geometry. The structural contribution to 
the Stokes shift is about 0.5 eV. The decomposition of the emission spectra in terms of independent- 
particle transitions is similar to what is found for absorption. For the emission, very weak transitions 
are found in Ge-rich clusters well below the strong absorption onset. 



PACS numbers: 71.15.Mb, 78.67.-n, 78.67.Bf 
I. INTRODUCTION 

The semiconductors silicon and germanium can form 
a substitutional solid solution of the form Gei-^Sia; cov- 
ering the whole range of compositions x. Pure Si is 
the most widely used material for electronic applications 
since many years and its fabrication technology is highly 
developed. However, the indirect band gap of bulk Si 
presents a problem for light-emitting applications. A so- 
lution that has been proposed to circumvent this prob- 
lem is nanostructurization of Si in structures comprising 
porous Sifi nanowires^ as well as Si nanocrystals. More- 
over, ample use has been made of the fact that Ge can 
easily be combined with Si in heterostructures. In addi- 
tion, Ge nanocrystals in a matrix of Si02, SiC, sapphire, 
or Si have been investigated by many groups^^ In many 
of those cases, intermixing between Ge and Si is found in 
the nanostructures. 

The two materials show different properties upon 
nanostructurization. While Si retains its character of an 
indirect material jSiLS the work of several groups showed 
that, depending on the structure, strain, etc., Ge nanos- 
tructures can become quasi-direct, i.e., they exhibit very 
strong transitions at or very close to the HOMO-LUMO 
transitioni&^&^ i 10 ! 11 : 12 The different behavior reflects the 
different character of the band gaps in the bulk materi- 
als. While both are indirect, the minimum gap in Si lies 
between T and a point near the X point, and the direct 
gap is much larger. In Ge, the minimum gap between 
r and L is energetically very close to the direct gap at 
r. The effects of confinement and structural relaxation 
result in a strong contribution of the T-T transition to 



the HOMO-LUMO transition, thus resulting in short ra- 
diative lifetimes of Ge clustersi&i£ 

Therefore, when mixing both materials two questions 
arise: what is the effect of the intermixing on the elec- 
tronic properties of the nanocrystals, and how are the dif- 
ferent characters of the two materials combined. While 
the answer to these questions is essentially well known for 
the bulk alloy, there are still few investigations concern- 
ing the mixed nanostructures. Most experimental studies 
use Stranski-Krastanov growth, (see, e.g., Ref . [lih . which 
results in relatively large structures in which it is safe to 
assume that the effects of confinement and alloying act 
independently However, for smaller structures this can- 
not be taken for granted. For small nanocrystals, pho- 
toluminescence experiments^ have been compared with 
theoretical results^ Furthermore, Ge-Si nanowires have 
been investigated experimentally^ and theoretically^ 

Previous theoretical studies have focused on the 
interplay of confinement and alloying i 16 i 18 However, 
these ab initio calculations were performed within the 
independent-particle approximation based on the Kohn- 
Sham scheme of static density- functional theory (DFT). 
Therefore, important many-body effects have been ne- 
glected, viz., the self-energy corrections describing the 
effect of the excitation of the electrons or holes individ- 
ually, as well as the electron-hole interaction. (Nonethe- 
less, these two effects are found to cancel each other to 
a large extent for many systems^) Furthermore, these 
calculations miss the very important depolarization or 
crystal local-field effects. 

The way we choose to improve upon the independent- 
particle results is provided by time-dependent 
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DFU22i2ii22i (TDDFT) in which the many-body ef- 
fects neglected in static DFT are introduced by the 
so-called exchange-correlation kernel f xc . With respect 
to the independent-particle results, the excitation 
energies are corrected, and the transitions between 
the independent-particle states are mixed. The de- 
gree of the mixing indicates the degree to which the 
independent-particle approximation fails to provide a 
good description of the system. 

Within TDDFT, we use the adiabatic local-density 
approximation (ALDA, also known as TDLDA) of the 
exchange-correlation kernel. The choice is motivated by 
the fact that TDLDA yields good results for optical spec- 
tra of isolated systems^ 3 - as well as for non-zero momen- 
tum transfers.— Note, however, that the ALDA is well 
known to fail in some cases, the most important of which 
are perhaps extended systems . 22 ' 25 Within TDDFT, the 
random-phase approximation (RPA) is obtained by ne- 
glecting the exchange-correlation kernel, i.e., by setting 
fxc = 0. From this, the independent-particle approxima- 
tion results from the neglect of the microscopic terms of 
the variation of the Hartree potential. In other words, 
the difference between the independent-particle approxi- 
mation and the RPA are the depolarization effects which 
are due to the inhomogeneity of the system. 

The simplest way to account for depolarization effects 
within an approximated classical picture is to apply the 
Maxwell-Garnett effective-medium theory to the complex 
dielectric function of the bulk alloy crystals. We applied 
the effective medium theory to our alloy clusters in or- 
der to check if the depolarization effects are dominant 
with respect to confinement effects and further many- 
body corrections. In fact, if it is the case, the spectra 
obtained using the effective-medium theory are in good 
agreement with full TDLDA results. These model cal- 
culations thus give a good overall description of the op- 
tical spectra at a much lower computational cost than 
TDLDA. 

In the present paper we first present TDLDA results 
for Gei-^Siz nanocrsytals, focusing on alloying effects. 

The results are then compared with the previous 
independent-particle calculations of Ref. [H, highlighting 
the effects of the depolarization, as well as of the mix- 
ing of transitions. We will also consider depolarization 
effects alone, decoupled from confinement effects, by ap- 
plying Maxwell-Garnett classical effective-medium the- 
ory. Finally, the emission properties of the nanocrystals 
are investigated by considering the geometries obtained 
after excitation of an electron-hole pair. 



II. CALCULATION DETAILS 

We used the same modeling scheme as used in Ref. [ID 
considering Gei-^Si^ nanocrystals with a fixed size made 
of 83 Gc and Si atoms. Quasi-spherical nanocrystals were 
built starting from one atom and adding nearest neigh- 
bors shell by shell, assuming bulk-like tetrahcdral coor- 




FIG. 1: (Color online) Example of the studied nanocrystal 
structures: Ge48Si3sHio8. The colors are: Si (yellow), Ge 
(blue), and H (white). 

dination. The outer bonds were saturated by hydrogen 
atoms. Alloying between Ge and Si was introduced by 
randomly exchanging Ge atoms by Si. The surface was 
then passivated with H atoms to saturate the remaining 
dangling bonds. 

In Ref. [T5, the study of Gei^S^ nanocrystals with the 
same number of atoms for ten different atomic configu- 
rations demonstrated that those with nearly uniformly 
distributed Ge and Si atoms possess the lowest total 
energies and nearly equal excitation energies. On the 
other hand, nanocrystals with deliberately clustered Si 
atoms and, hence, rather different excitation energies 
give rise to total energies substantially higher than the 
average. As their probability of occurrence is, conse- 
quently, small, the configurational average can be re- 
placed by the study of only one nanocrystal with nearly 
uniformly distributed Si and Ge atoms for each composi- 
tion x. We selected nanocrystals with a total number of 
83 atoms of Ge and Si, having a diameter of about 1.5 nm. 
This radius corresponds to a sphere of the volume occu- 
pied by 83 atoms in the bulk. These nanocrystals are 
large enough to exhibit the characteristics of nanocrys- 
tals as opposed to much smaller structures which show a 
molecular-like behavior i 6 i 13 i 26 Moreover, they are small 
enough to present significant confinement effects on the 
electronic states. An example of the atomic arrangement 
is shown in Fig. [T] for Ge4sSi35Hi 8. 

To obtain the relaxed geometries of the Gei-zSLj 
nanocrystals, we used the plane-wave code VASP— 
within the local-density approximation (LDA) in the 
parametrization of Perdew and Zungcr 2 * and the 
projector-augmented wave (PAW) method^. This com- 
putational set-up is the same as the one employed in 
RefG3 

Starting from the relaxed geometries we obtained 
the optical spectra at zero temperature using TDDFT 
as implemented in the computer code octopusi 30 i 31 
The electron-ion interaction is described through norm- 
conserving pseudopotentials 3 - 2 - and the LDA 2 ^ is em- 
ployed in the adiabatic approximation for the xc poten- 
tial (i.e., we applied the TDLDA). The time-dependent 
Kohn-Sham equations, in this code, are represented in 
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FIG. 2: (Color online) Absorption spectra as a function of the Si content of the clusters: Independent-particle spectra (dashed 
lines) compared with TDLDA results (solid lines). The green arrows (arrows at lowest energy) mark the HOMO-LUMO gap, 
while the red arrows (middle arrows) mark the ASCF excitation energies and the blue arrows (arrows at highest energy) mark 
the first transitions of the Casida's analysis. The independent-particle curves are divided by a factor of 15. Note the different 
scale in the two panels. 



a real-space regular grid, using a spacing of 0.275 A at 
which the calculations are converged. The simulation box 
is constructed by joining spheres of radius 4.5 A, centered 
around each atom. 

To calculate the optical response we excite the system 
from its ground state by applying a delta electric field 
EoS(t). The real-time response to this perturbation is 
Fourier transformed to get the dynamical polarizability 
a{u>) in the frequency range of interest. The absorption 
cross section cr(u>), is then obtained from the relation: 

4-7T 

cr(u)) = — u) 3m auo) , (1) 
c 

where c is the velocity of light in vacuum. For a tech- 
nical description of this method we refer to Refs. 30 
and [H. To obtain the RPA spectra in this approach we 
just kept the exchange-correlation potential fixed during 
the propagation (which amounts to making / xc = 0). 
The independent-particle spectra were calculated by fur- 
ther fixing the Hartree potential to its initial value. A 
time step of 0.0055 ft/cV and a total propagation time of 
37.5ft/cV were sufficient to ensure a stable propagation. 
We estimate our numerical precision in the spectra to be 
better than 0.05 eV. The mixing of independent-particle 
transitions in the spectra has subsequently been investi- 
gated using the Casida's formulation of linear-response 
TDDFT.Si^ 

In order to test the coherence of our calculations, we 
compared independent-particle spectra computed using 
matrix elements calculated with the VASP code£ ' 36 i 37 to 



spectra obtained by time propagation by means of the 
octopus code. No substantial differences were found be- 
tween the results. 

Due to the interest in the luminescence of this kind 
of systems, the lowest excitation energies are of partic- 
ular importance. As a consequence, we decided to com- 
pare three different quantities. Besides the lowest ex- 
citations of the TDLDA absorption spectra, we present 
the Kohn-Sham HOMO-LUMO gap of the ground state 
and the ASCF excitation energies calculated as the dif- 
ference Sascf = E(N,e-h) - E(N), where E(N) is the 
total energy of the ground-state and E(N, e-h) is the to- 
tal energy of the excited-state configuration where one 
electron has been promoted from the HOMO into the 
LUMO. In this way, an excited configuration is mod- 
eled and the electron-hole interaction is partially taken 
into account. This approach enables also a subsequent 
ionic relaxation with the electron-hole pair present, which 
yields the description of the excited-state geometries. Us- 
ing the excited-state geometries we were also able to eval- 
uate emission spectra. 



III. RESULTS 

A. Ground-state geometries: absorption 

Results for the photoabsorption of the Gei-^Si^ 
nanocrystallites are presented in Fig. [21 In the figure we 
compare the independent-particle response (dashed lines) 
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with full TDLDA calculations (solid lines) for the whole 
range of compositions. The dependence of the curves on 
the composition turns out to be quite smooth. When go- 
ing from Ge-rich clusters to Si-rich clusters, we observe a 
shift to higher energies of the onset and a suppression of 
the absorption strength of the first peak (note the differ- 
ent scales of both panels). Moreover, the dependence 
of the peak position on the composition x is roughly 
the same in the independent-particle and the TDLDA 
schemes. For the clusters studied, the total intensity of 
the independent-particle spectra is strongly suppressed 
in TDLDA (the independent particle curves are divided 
by a factor of 15 in Fig. [2]). Indeed, this quenching of the 
absorption is a well known effect that is due to the inclu- 
sion of classical depolarization effects, and not from the 
exchange-correlation effects accounted for by the TDLDA 
kernel. This can be verified by calculating the absorption 
spectrum within the RPA, which includes all classical ef- 
fects due to the variations of the Hartree potential but 
neglects quasi-particle and excitonic effects. Indeed, the 
spectra wc calculated within RPA are so close to the 
TDLDA spectra shown in Fig. [2] that we chose not to 
show them. 



B. Transition analysis and excitation energies 

For Ge-rich nanoparticles, both the position and the 
composition dependence of the peaks are already well 
described at the level of the independent-particle approx- 
imation. This perhaps surprising fact can be explained 
in terms of compensation of quasiparticle corrections and 
binding energies of the excitons. For small x at the ab- 
sorption edge, the first peak is strong and appears essen- 
tially at the HOMO-LUMO energy gap. The absorption 
edge has a completely different nature in Si- rich clusters. 
In fact, in this case the peaks of lowest energy have a 
vanishing oscillator strength in TDLDA, thereby blue- 
shifting the absorption edge with respect to the HOMO- 
LUMO gap. 

In order to analyze the origin of the different peaks 
in the spectra as a function of the composition of the 
Gei-^Sia, alloy, we decomposed the excitations in sums 
of Kohn-Sham particle-hole transitions through the so- 
lution of Casida's equation . 34 ' 35 We found that on the 
Ge side of the composition range, the lowest peak of 
the spectra which defines the absorption edge is pro- 
duced essentially by a strong, pure transition between 
the Kohn-Sham HOMO and LUMO. This is certainly one 
more reason for the similarity between the independent- 
particle and the TDLDA spectra. The large peak with 
HOMO-LUMO character decreases in intensity with in- 
creasing percentage of Si in the Gci-^Si^ alloy, until the 
composition of about x = 0.2 when it disappears. For 
even smaller x the absorption at the onset is determined 
by a strong mixture of transitions between states close 
to HOMO and LUMO. However, the very lowest transi- 
tions are forbidden, producing a significant blue shift of 
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FIG. 3: (Color online) Absorption spectra as a function of the 
Si content of the clusters: Effective-medium theory (dashed 
lines) using the dielectric function of bulk Gei-^Si^ , com- 
pared with TDLDA results (solid lines). 

the absorption edge. In the intermediate energy range 
and for all compositions, excitations can be decomposed 
as a sum of many contributions. 

In Fig. H] we also show the HOMO-LUMO gap (green 
arrows), the lowest excitation calculated within the 
ASCF approximation (red arrows), and the first exci- 
tation within TDLDA (blue arrows, these excitations are 
dark on the Si side of the composition). For all cases 
we find that the first ASCF excitation is, as expected, 
blue-shifted with respect to the HOMO-LUMO gap, and 
that the first TDLDA transition is at slightly higher en- 
ergies. The differences are, however, quite small, and of 
the order of one tenth of an eV. This is known for this 
class of systems^ part of it being due to the cancella- 
tion of self-energy and excitonic effects^ However, the 
differences appear to be slightly larger in the region of 
intermediate composition x, i.e., of a greater degree of 
structural disorder. Similarly, Degoli et al. found in Si 
nanocrystals that the differences become larger in cases 
of stronger localization^ It is also clear from the plot 
how the first transition becomes forbidden while going 
from Ge-rich to Si-rich clusters — this behavior is remi- 
niscent of the different character of the band gap in the 
parent bulk Ge and Si. 



C. Absorption from classical effective- medium 
theory 

A much simpler approach to model the absorption 
cross section of a Gei_ x Si s nanocrystal is to start from 
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the complex dielectric function of the corresponding bulk 
alloy crystal and to apply the effective-medium the- 
ory.—^ This classical approach is based on Maxwell's 
equations and neglects completely the microscopic de- 
tails, such as atoms and bonds. Of course, this assump- 
tion is better justified when the size of the system is large. 
However, it always handles correctly the boundary condi- 
tions for the Maxwell's equations at the interfaces, which 
give the very important contributions to the dielectric 
response through the classical depolarization effects. Of- 
ten, these classical contributions are enough to describe 
the physics of the dielectric response of a composite sys- 
tem made of objects embedded in some matrix i 41 ' 42 i 43 
Our clusters can be considered as a family of spheres 
of volume Vobj cut from a Gei-^Si^ bulk alloy. The 
Maxwell-Garnett expression^ in the specific case of an 
isolated spherical object in vacuum yields: 



/ * W T/ 3m e(o;) 

a(u!) = 9 — l/ bi n 

c J [<HeeH + 2] 2 + [ 



3m e(u>)Y 



(2) 



where e is the experimental complex dielectric function 
of the bulk alloy. To represent well the extension of the 
polarizable nanocrystals, we took an average distance of 
the furthermost saturating hydrogen atoms to obtain the 
radius of the cluster, and consequently Vobj. This results 
in a radius slightly larger than the one mentioned above 
which takes into account only the Ge and Si atoms. We 
used the value of i? Q bj = 9.1 A, but the results are fairly 
insensitive to a (reasonable) choice of this value. The ex- 
perimental dielectric function of the bulk alloy with pre- 
cisely the needed composition x has been obtained from a 
discrete set of measurements^ using a recently developed 
interpolation schemed This scheme interpolates 3m e(uj) 
making use of the screening sum rule and the positive 
definitcness of the spectra. £Re e(w) was subsequently ob- 
tained by means of the Kramers-Kronig relations after 
fitting appropriate tails to the imaginary parts. The lat- 
ter procedure was tested on the input curves to insure 
the quality of the real part. 

In Fig.[3]we show spectra calculated using the effective- 
medium theory for selected compositions. Comparing to 
the first-principles TDLDA curves, we can see that, al- 
ready for this relatively small size of clusters, the classical 
theory gives a quite good overall description of the ab- 
sorption spectrum. As expected, effective-medium the- 
ory is not capable of describing the peaks of the individ- 
ual transitions, but it describes correctly the intensity 
and the trends of the spectrum. Once again these re- 
sults confirm that the dependence on the composition of 
the optical spectra is smooth and that the confinement 
and alloying effects act independent to a large degree. In 
fact, within this classical scheme, only the alloying ef- 
fects determine the variation of the absorption response 
as a function of x. In particular, the confinement-induced 
opening of the HOMO-LUMO gaps and the resulting 
blue-shift of the absorption onset are not accounted for. 
The confinement effects could be described by introduc- 
ing a size-dependent crystallite dielectric function which 



can then be used to calculate the spectra of the crys- 
tallites in a different environment^ In this sense, com- 
parison of the TDLDA results with the present effective- 
medium approach gives an idea of the importance of the 
confinement effects on the overall spectra. 

We note that these model calculations can be per- 
formed at negligible computational cost, and therefore 
provide a simple and fast method to obtain reasonable 
spectra for medium and large nanocrystallitcs. 

Given the strong influence of the depolarization effects, 
the question arises as to why the independent-particle 
spectra have been successfully compared previously with 
experiment^ This agreement is, in fact, not a fortuitous 
coincidence, but due to the experimental conditions. The 
experiment has been done on Ge NCs inside a matrix 
of sapphire4£ This reduces strongly the depolarization 
effects because it reduces the inhomogeneity of the sys- 
tem. The calculations, on the other hand, treated NCs in 
vacuum but neglected the depolarization effects. There- 
fore, together with the cancellation between self-energy 
effects and electron-hole interaction mentioned above, 
the independent-particle approximation provides in fact 
a good description of the spectra of this particular exper- 
iment. 



D. Excited-state geometries: emission 

In order to calculate the emission properties, we con- 
sider the geometry of the relaxed nanocrystals where 
the ionic relaxation has been carried out after transfer- 
ring an electron from the HOMO to the LUMO Kohn- 
Sham orbital. As the radiative lifetimes are usually much 
longer than the times that the electrons (holes) take to 
relax to the LUMO (HOMO), we can assume thermal- 
ized electron-hole pairs. As their lifetimes are then de- 
termined by the exponential factor describing their dis- 
tribution^^ it is the onset of the emission which re- 
flects the emission properties, while the higher parts of 
the spectrum are suppressed. 

Stimulated emission spectra can be easily obtained 
within our formalism by calculating the absorption cross- 
section (T a bs(w) at the excited-state geometry. Lumi- 
nescence spectra can then be calculated from the van 
Roosbroeck-Shockley model 4 ^ 
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?hcj/k h T ^ 



5"abs(^) , 



(3) 



where fcb is the Boltzmann constant and T is the tem- 
perature. 

In Fig. [4] we compare the absorption cross sections 
for the ground-state (solid lines) and the excited-state 
(dashed lines) geometries for Gei-zSLj clusters in all 
the composition range. The energy difference between 
the onsets of absorption and emission corresponds to the 
structural contribution to the Stokes shift. 

In Fig. O we compare the independent-particle re- 
sult with the TDLDA for the absorption <7 a bs(w) at the 
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excited-state geometry. The conclusions with respect 
to their similarity drawn for the ground state remain 
valid for the excited-state geometries. The same ap- 
plies to the transition analysis using Casida's equation. 
The character of the emission onset on the Ge side of 
the composition range is already well described within 
the independent-particle approximation. The decompo- 
sition of the excitations as a sum of Kohn-Sham tran- 
sitions provides a picture strictly analogous to the one 
for the absorption spectra: the lowest transitions of the 
Ge-rich nanocrystals correspond to almost pure Kohn- 
Sham transitions, while for large x and for the higher 
transitions, independently of x, strong mixing is found. 

However, it is important to note that a very weak peak 
appears at about 1.95 eV. This is red-shifted by about 
0.5 eV with respect to the ground-state calculation. The 
peak, which can be easily seen in Fig. [5J is clearly present 
in both the independent-particle result and the TDLDA 
result. It occurs for all compositions on the Ge side of the 
x range and is almost composition independent. It ap- 
pears to be connected with the lowering of the symmetry 
as compared to the ground state where the pure Si or Ge 
nanocrystals without alloying have T^ symmetry and the 
HOMO-LUMO transition is threefold degenerate. The 
symmetry breaking due to the geometry relaxation un- 
der excitation has therefore a much stronger effect than 
the introduction of the alloying, which at the Ge-rich 
side splits the degeneracy only slightly and which does 
not change the character of the strong absorption onset 
at the HOMO-LUMO transition. Due to the argument 
above, cf. Eq. [3J the appearance of this weak peak will 
strongly increase the radiative lifetimes of the systems. 
We conjecture that this might be responsible for the fact 



that even though several theoretical predictions coincide 
in that Ge nanostructures should have strong transitions 
at the absorption onset, few experiments have been able 
to detect luminescence from excitons in Ge nanocrystals. 



IV. CONCLUSIONS 

The absorption spectra of free Gei-^Si^ have been 
calculated within time-dependent density-functional the- 
ory in the adiabatic local-density approximation. The 
changes of the spectra upon changing composition x are 
smooth. In particular at the absorption onset, the po- 
sition and the composition dependence of the spectra 
is found to be already well represented by independent- 
particle results. The analysis of the solutions of Casida's 
equation show that this is due to the fact that the first 
transition of Ge-rich nanocrystals corresponds to an al- 
most pure independent-particle transition. The TDLDA 
onsets are slightly blue-shifted with respect to their 
independent-particle counterpart, their composition de- 
pendence at the Ge side of the compositional range is 
practically the same. For higher Si contents, a mixing of 
many independent-particle transitions is found. 

Depolarization effects are strong and their inclusion 
alone, even within a simplified classical model, on top of 
an independent-particle calculation allows to get the cor- 
rect physical picture of the optical response. They can 
be approximately taken into account at a negligible com- 
putational cost using Maxwell-Garnctt effective-medium 
theory. 

Further many-body terms do not modify significantly 
the spectra due to the cancellation of opposite contri- 
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FIG. 5: (Color online) Absorption (red) and emission (black) spectra of the pure Ge nanocrystal in the independent-particle 
approximation (dashed) and the TDLDA (solid). Note the different scales. 



butions given by quasiparticle corrections and excitonic 
effects. As a consequence, the effects of the excitation 
of an electron-hole pair do not alter the comparison be- 
tween the independent-particle spectra and the TDLDA 
results. 

Emission spectra have been investigated using the ge- 
ometry of excited nanocrystals. A Stokes shift of about 
0.5 eV is found. Very weak peaks appear at the absorp- 
tion onset for all systems on the Ge side which suppress 
the radiative transition probability and lead to long ra- 
diative lifetimes. 
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